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Abstract. The large multiplicities at the LHC may permit flow harmonics to be determined 
on an event by event basis in Pb-Pb collisions. We extract these harmonics from event by 
event di-hadron correlations. Within a fine centrality bin, we find the correlation function varies 
substantially on an event by event basis, indicating large fiuctuations in the initial conditions for 
a given impact parameter. Such large fluctuations lead to some events being highly triangular 
or highly elliptical, where the angular correlation function is completely dominated by the 
respective second and third Fourier harmonics. We will show unfolded V2 distributions for 
various centralities, and implications for our understanding of the initial conditions. 



1. Introduction 

Relativistic heavy-ion collisions aim to create the QGP (Quark-Gluon Plasma), a unique state of 
matter where quarks and gluons can move freely over large volumes in comparison to the typical 
size of a hadron. Spatial inhomogeneities in the initial state can create pressure gradients in the 
QGP, which induce angular correlations known as azimuthal flow. It is convenient to express 
the single particle density of produced particles as a Fourier series [I]: 

— OC 1 2 ^ VnCOs{n{<p - Tpn)) (1) 
^ n=l 

where Vn and ipn are the magnitude and direction of order n harmonic flow. The magnitude 
Vn is typically extracted from azimuthal correlations averaged over an event ensemble, where 
the second order (elliptic flow) is usually the dominant contribution. It has been proposed that 
event by event fluctuations in the initial state give rise to the odd orders [2]. They also lead to 
variance in Vn (ctjn) for all orders, that can be extracted from data via two and four particle 
cumulants [B]. These fluctuations should be also evident in the event-wise pair distributions. 
All orders of Vn can be extracted without the need to identify ipn using the pair density: 
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where A(j) is the angular difference between the two particles in a pair. For factorization, it 
is important that the only correlation between the two particles is due to flow. In this case, 
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Figure 1. Left Panel: Event averaged Ac/) distributions for Pb-Pb y/sJjN = 2.76 TeV collisions. 
The dot-dash lines show the sum of a Fourier decomposition (with n = 2 and n = 3 contributions 
below). Right Panel: Event wise VnA extracted from 50 events in chronological order. Dashed 
lines show ztlo" statistical fluctuations for events with true w„ = 0. 



VnA factorizes into the product of Vn, VnA = 'v'^- Factorization has been demonstrated for 
these energies at low [3\. Measurements of di-hadron correlations from ALICE show V2 
and f3 manifest themselves as a nearside "ridge" and awayside "double hump" for very central 
collisions O 01 [5] . We present measurements of event by event di-hadron correlations in Pb-Pb 
y^ijviv = 2.76 TeV collisions. The large multiplicities at the LHC and near full pT acceptance 
of the ALICE detector favour such measurements. Values of V2A and V3A are extracted from 
these distributions, and we will show events with large values of V2A and V^a- Measurements 
of event wise V2A allow for the full V2 distributions to be obtained, which will also be shown. 



2. The correlation function 

We use charged hadrons with transverse momentum px > 0.15 GeV/c from the ALICE TPC to 
form the pair correlation function. This is defined as follows: 

C(A0), = i^^EilMi (3) 

\-^^ pairs) 

where Npairs{A(j))i is the number of pairs in bin i, and {Npairs) is the average number of pairs 
summed over all bins. It therefore quotes the relative deviation from the mean number of pairs. 
The main advantage of the correlation function is that pair inefficiencies due to the detector (e.g. 
tracking inefficiencies) in the numerator and denominator cancel to leading order. To extract 
VnA from C(A(/)), a definition used in a previous ALICE publication will be adopted [3]. It is 
given by: 

VnA = E CiCosinA<p)/ ^ Ci (4) 



where n is the harmonic of interest and Ci is the value of the correlation function in a particular 
A(p bin. The key advantage of this method is that the calculation of VnA is invariant with 
respect to the choice of normalisation of the correlation function. The left panel of figure [1] 
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Figure 2. Ac/) distributions from events 14 and 46 in the right panel of figure [TJ 



shows the event averaged correlation function obtained for Pb-Pb centrality 4-5%. A Fourier 
decomposition is performed, and each coefficient is extracted using equation [H The Fourier 
series is shown by the black hnes. It was found that the first 5 terms are needed for a good 
description of the data for this pT range, and one can obverse the n = 2 contribution is larger 
than the n = 3 contribution. 



3. Event by event Vat a 

Rather than obtaining C{A(f)) and VnA from an event ensemble, we investigate these 
measurements in single events. In the right panel of figure [H event-wise values of V2A and V3A 
are shown for 50 events that occur for centrality of 4-5% chronologically. The dashed lines show 
expected itlcj levels of background fiuctuations for each order, due to statistical fluctuations. 
There are clearly many events with V2A higher than expected from just background fluctuations. 
Such an observation indicates large fiuctuations in the initial conditions, even within a small 
centrality range. The V2A fluctuations are bigger than the V3A fluctuations, and this is consistent 
with the fact (F2A) is bigger than (Vsa) (shown in the left panel). This is also inconsistent with 
background fluctuations, where we find all orders of n fiuctuate equally (the red and blue dashed 
lines overlap). In figure El the A(p distributions for event 14 and event 46 from the left panel 
of figure [1] are shown. It is observed that these events are completed dominated by elliptical 
(n = 2) and triangular (n = 3) flow respectively, and are in contrast with the distribution seen 
in the left panel of flgure [TJ 



4. Unfolded V2 distributions 

In flgure El we show V2 distributions as a function of centrality. These are obtaining by unfolding 
raw V2A distributions and applying the transformation /(f2) = fi\/^2A), where /(1^2a) is the 
unfolded V2A distribution. The unfolding procedure is needed to remove the effects of the 
non-flow and trivial statistical smearing. The statistical smearing is determined by sampling a 
dN/dcj) distribution with a known input V2A) then extracting the event- wise reconstructed V2A- 
The sampling is repeated, so one obtains the distribution of reconstructed V2A for the given 
input. The contribution of non-flow is determined by fitting the pseudo-rapidly dependence of 
V2A with a Gaussian (non-fiow) and constant (fiow). The systematic uncertainties (shown by 



Pb-Pb ifs^ = 2.76 TeV 



10 I 

: Method: = V^^, > 0.15 GeV/c 

' T*>- - - Bessel-Gaussian fit 

1 n \ 



10 



lOV 



10- 



1 • 



,11 



1 



RLICE 

PRElIMI^flRY 



*0-5% 
♦5-10% 
* 10-20% 



10 



20 



30 

V2 (%) 



10' 



110" 



10' 



10 



1 - 



10-^ 



lO" 



* 20-30% 
♦30-40% 
* 40-50% 




x-^... ALICE 

PHELIMINflRY 



10 



\ XM 



^a** - 



20 



30 

V2 (%) 



Figure 3. Unfolded V2 distributions as a function of centrality, subject to a Aij > 0.8 gap. 



the boxes) are due to the uncertainty in the non-flow contribution (evaluated by changing the 
Atj cut), statistical smearing hypothesis (evaluated by changing the number of particles used) 
and centrality selection. The V2 distributions are expected to reflect the participant eccentricity 
distributions [8], therefore provide constraints on models of the initial conditions. We fit these 
distributions with a Bessel-Gaussian distribution, which is the relevant p.d.f if the components 
of the eccentricity vector {ex,Sy) are subject to Gaussian fluctuations [9]. These fits generally 
describe the data well, bar central collisions. We found (not shown) that the unfolded (^2) and 

\/ (^2) + 2*^^ consistent with ALICE measurements of f2{4}, and V2{2} respectively [4J, which 
is expected if the f 2 distributions are Bessel-Gaussian. The widths of the unfolded distributions 
were also found to be consistent with the widths obtained via ^/0^5(v2{2y^'^^^J2{^y^ ■ 

5. Summary 

We have investigated event by event di-hadron correlations, and found some events dominated by 
n = 2 (elliptical events), and others dominated by n = 3 (triangular events). The observations 
are beyond what is expected from trivial statistical fluctuations. We have presented unfolded V2 
distributions and found they are generally described by a Bessel-Gaussian distribution. Theses 
distributions can be used as constraints for models of initial state eccentricities. 
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